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We compare recent CEBAF data on inclusive electron scattering of 4.05 
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GeV electrons on nuclei with predictions, based on a relation between struc- 

0^ 



ture functions (SF) of a nucleus, a nucleon and a nucleus of point-nucleons. 
The latter contains nuclear dynamics, e.g. binary collision contributions in 
^ ' addition to the asymptotic limit. The agreement with data is good, except 

c ■ 

in low-intensity regions. Computed ternary collision contributions appear too 
^ I small for an explanation. We perform scaling analyses in Gurvitz's scaling 



variable and found that for yc ^ 0, ratios of scaling functions for pairs of 
nuclei differ by less than 15-20% from 1. Scaling functions for {yc) < are, 
for increasing Q^, shown to approach a plateau from above. We observe only 
weak Q^-dependence in FSI, which in the relevant kinematic region is ascribed 
to the diffractive nature of A^A^ amplitudes appearing in FSL This renders 
it difficult to separate asymptotic from FSI parts and seriously hampers the 
extraction of n{p) from scaling analyses in a model- independent fashion. 
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In the following we discuss some aspects of inclusive scattering of high-energy electrons 
from nuclear targets with the following cross section per nucleon 



dVLdv A 



aMiE; 6, v) 



(1) 



Q 

The inclusive, as well as the Mott cross for point-nucleons aM, are measured as functions of 
the beam energy E, the scattering angle 9 and the energy loss z/. The above F^2{^j Q^) are 
two nuclear structure functions (SF) which account for the inclusive scattering of unpolarized 
electrons from randomly oriented targets. These depend on the squared 4-momentum = 
— and the Bjorken variable a;, corresponding to the nucleon mass M with range 
< Q'^jlMv < A. 

We concentrate on the recent CEBAF E=4.05 GeV inclusive scattering data on nuclear 
targets, notably Fe0. These cover a range of squared 4-momenta transfer = — ■, 1 < 
(52(GeV^) < 7 and Bjorken variable x = Q'^/2Mu, 0.20 < x < 4.2, and vastly extend the 
angular and energy-loss ranges of the older NA3 SLAci and related experimentsi@. 

Predictions for the above CEBAF experiment have been made beforeH. These were based 
)del which connects the structure functions (SF) of a nucleus and of a nucleoni 



on a moc 



F,\x^Q') = l (2) 

A calculation of the above nuclear SF requires two input elements. The first one is 
F^ = F^^'"^, properly weighted with proton and neutron numbers in the nuclear tar- 
get. The F'^ contain nucleon elastic (NE) part in terms of the standard static electric 
and magnetic form factors G^, G^, and nucleon inelastic (NI) parts which are available in 
parametrized and tabulated formsB. The second element is the structure function f^^ for 
a nucleus of point-nucleons, which accounts for the nuclear dynamics in F^. Calculations 
employ a relativistic extension of the non-relativistic (NR) Gersch-Rodriguez-Smith (GRS) 
series in 1/q i. The latter contains an asymptotic limit (AL), related to the single-nucleon 
momentum distribution n{p), and Final State Interactions (FSI) which are dominated by 
binary collisions between the knocked-on nucleon and a nucleon from the core. 



Expressions, similar to (|]), have previously been suggested with light-cone momentum 
fractions instead of Bjorken variablesi'0. The two coincide for Q"^ ^ oo and Eq. 
conjectures the approximate equality to hold for Q^, in excess of some Q^. 

The model incorporated in (2) locates weak A-dependence of Fj^{x, Q"^) only in the 
neutron excess SN/2A, and in f^^ i, thus 

F,^{x,Q')^F,ix,Q') + Oil/A); A>12, (3) 

with the above restriction on A, due to the neglect of nucleon recoil in calculations of 
in (0). In Fig. 1 we compare Fe dat J and the above predictions! and conclude: 

i) There is good agreement in the deep-inelastic region u > Q'^/2M, and satisfactory 
correspondence on the immediate NE side of the quasi-elastic peak (QEP), x > 1. 

ii) Data are for given 6, and not for fixed Q^, but from the fact that increases with 
6, and for given 6 with decreasing z/, one observes that discrepancies grow with decreasing 
6, i.e. for decreasing Q^. For < 1.6 GeV^ this may in part be due to variation in some 
ill-determined A^A^ scattering parameters, which enter the FSI calculations. However, it is 
more likely, that < Ql{x,6), below which the representation becomes progressively 
flawed. 

iii) For all 6, computed cross sections overestimate the data by a factor up to 2-3 for the 
lowest energy losses u, where cross sections have dropped orders of magnitude from their 
maximum. 

In an attempt to understand the above discrepancy we worked out more complex FSI 
contributions, generated by ternary collisions (TC). These are on general grounds expected 
to effect low-intensity regions far from the QEP, where they compete with both the AL and 
FSI, due to BC. Their small, but noticible effect has recently been established in the response 
of liquid ^He 0. We have used a relativistic extension of the above TC contribution and 
found that these indeed contribute in the above kinematic regions, but only with insignificant 
weight. 

At this point we recall that it is not at all clear that the above is a real discrepancy. Using 



different acceptable n{p) produce results which range over the area of the above mentioned 
local discrepancies (see Figs. 5,6 in Ref. |^). 

The NA3 experimenti has also been analysed by means of versions of the Plane Wave 
Impulse Approximation (PWIA) in terms of a spectral function, occasionally supplemented 
by additional FSlEHll. We recall attempts by Benhar et al to ascribe to color transparency 
a desired lowering of cross sections^. Conventional 2p — Ih FSI on the PWIA apparently 
produce the same effect^. 

A comparison with the above mentioned generalized GRS results shows excellent agree- 
ment in inelastic regions. For decreasing energy losses Ciofi and Simula somewhat underesti- 
mate intensities0, whereas our approach overestimates them. In fact, one wonders how two 
entirely different series, GRS and lA tend to produce comparable answers. The explanation 
is a recent demonstration for NR dynamics of a surprising, order-by-order correspondence 
in l/q0. 

Next we turn to a few selected scaling analyses, previously applied to the NA3 data 
We had used there a relativistic West-GRS scaling variable suggested by GurvitJll yc, 
related to x by ?/g ~ {Mh'/q){l — x). We investigated ratios of inclusive cross sections 



and in particular ^^'^^^{yc < in the NE region. For the latter, data for selected 
nuclei, nuclear matter and in fact for all investigated targets treated simultaneously, showed 
universal coarse scaling in yc- The above then also holds for the ratios ^^"^'^^{yc < 0, Q"^). 
The same resultsc3 when a relativistic Fermi gas scaling variable is used§. 

We report below on such an analysis of the new CEBAF data for the pairs C,Fe and 
Fe,Au. A novel element is the extension of the yc-Tange into the entire NI range yc > 0. 
Although for all 6 or Q^, cross sections per nucleon for given yc span 4-5 decades (or about 
3 decades for ^^'<^>)^ ^Ai,A2 ^ within 15-20% and frequently better (see Table I). This 
agrees with the approximate A-independence (0) of nuclear SF. Occasional larger deviations 
from 1 are ascribed to experimental uncertainties, in particular for data of lowest intensity. 



Next we consider the EMC ratio ^ which is not a special case of ^'^^'^'^ with A2 — »• 
(A^) = D/2. Density, momentum distribution and pair-distribution function for D differ and 
A > 12. Although the NE and the NI regions both contain information on the single-nucleon 
momentum distribution n{p), implicit in f^^, Eq. (|^), simple expressions for ^^'^ can only 
be given for yc < 0. However, that portion is not directly observable. In particular for high 
Q^, NI contributions compete with the NE ones, even on the elastic side 1/g ^ (x > 1) 
of the QEP (cf. Ref. |^, Fig. 4), as has already been realized in the analsis of early high-E' 
experiments on the lightest nuclei D, ^ He, ^Heil. 

In order to isolate the NE part of (."^''^{yc < 0, Q"^), one has to remove NI parts from the 
data, which can be done in several ways. For instance (|^) in conjunction with (1), provides 
a model for the NI parts, which gives perfect agreement with the data for ?/g ^ 0- Although 
this need not be the same for < 0, we assume this to be the case. The procedure 
ultimately becomes impractical, because both total and NI parts rapidly decrease with ya 
beyond yc ~ —0.25. One has then to rely on directly calculated NE parts, using again (||), 
now with NE parts for . The result is 



Fr''\x,Q') = l[G%{Q')rS{x-l) 

Fn,Q^) - mOT^mnis^^ _ 1) (5a) 
Fr'\x,Q') = lr''ix,Q')[GUQ')r 



l + r] 

Substitution in (1) yields for Eq. m [t] = Q'^/AM^ 



(x2mVg^)([Gf ]2 + r/K]2)/(l + r/) + tan2(^/2)([G^]2i 



ei^'^ {x ^ 1, Q2) = fPN[x ^ 1, Q2) (6b) 

where for better readability, we have dropped the arguments on . 

In Fig. 2 we plot ^ne^ iVG < 0,(5^) against for a number of narrowly binned yc- 
Whenever possible, we give the two results, which with the exception of yc = 0, approx- 
imately agree. For all ?/g < in the kinematic range of the experiment, C,ne approaches 



the AL — s> oo from above, or can confidently be extrapolated to beyond the observed 
ones. For the largest there is hardly any Q"^ dependence. 

We mention here an analysis of a PWIA scaling function which differs from the above ^"^'^ 
primarily by a kinematic factor and by the use of a different scaling variable uq = y^^^^ 00. 
Arrington etal found that for i/q < —0.3 GeV, ^^^^^^ approaches a plateau from above, 
with persistent structures for large Q^. However, for yo = 0.0, —0.1, ^^^^^^ increases with 
growing Q^. Assuming that NI contributions have been properly removed, the spelled-out 
differences may be due to the use of different scaling variables. Those amount to the implicit 
retention of different FSl parts which are not easily compared^. Since the authors of Ref. 
|l] plan studies in still different variables, we do not pursue the issue here. 

A plateau is conventionally related to the AL, from which one extracts the single-nucleon 
momentum distribution n{p). However, the standard argument becomes invalid, if parts of 
the FSI happen to be weakly Q^-dependent, causing the plateau to also contain FSI parts. 
There are strong indications that this is approximately the case for the kinematic range of 
the CEBAF experimentii, which seems to contradict the behavior of ^^'^ for low Q^, seen 
in Fig. 1. 

We give here only an outline of the resolution of the apparent contradiction: details may 
be found in Ref. |^. The argument follows the construction of f^^{x, Q"^) from the NR SF 
4>{q,yG)- For smooth underlying A^A^ interactions the g-dependence of FSI may as follows 
be written and parametrized 

y) = Fo{y) + ^'^iq, y) + A--0(g, y) (7a) 
Fo(?/) = lim (p{q,y) = — / dppn{p,'jk), (7b) 
^"'\q,y) = U^"\q)yY.al:^y'''exp{-[A("^y]') (7c) 

n 

r^-{q,y) = t/(^)(g) E4^^Z/'"exp(-[A(^)y]2) ^^^^ 

n 

Fo and A0 are parts even in y, 0°'^°' is odd. 

The notion of a local, energy- independent Vat at breaks down for increasing lab momentum 
q, and V has to be replaced V Vejf = tq, the latter being the off-shell elastic T matrix 
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for small momentum transfer. It is the dominant lm[t{q, = 0)] oc which is hardly 

g-dependent in the relevant range. The link of with f^^ in (Q) rea dJ 



lim D(x,Q^) = l, (8) 



with D a purely kinematic factor. From (pEf) one concludes that for ?/ ^ {x > 1) 



^NE iVG < 0, Q') ^ f^ivG < 0, Q') = MD<P{q, VG < 0) (9) 

Consequently a weakly g-dependent 0(g, i/g) acquires dependence through D for relatively 
low Q^, which gets weaker for increasing Q, as indeed observed. 

It is of interest to recall here a suggested parametrization of C^^(2/o, Q"^), function of 

e^iyo) = + C,exp(-%o|) (10) 

tt^ +1/0 

The above appears to provide a universal fit for g-dependent data, with practically A- 
independent parameters for A > 12, which follows from Eq. (0). Moreover we note that an 
expansion of (0) in yo leads to a form similar to ([Td]). Without entering an interpretation 
of the parameters, this implies that either all the data are in the AL region, or that FSI in 
addition to the AL are hardly g-dependent. If the latter is the case, the extraction of n{p) 
from the AL, without some knowledge of dynamics, becomes close to impossible. 

Summarizing we put into focus the relation (0) which enables the determination of nuclear 
SF for computed SF f^^ of a nucleus, composed of point nucleons. The latter contains 
the AL, expressed by the single-nucleon momentum distribution n{p), and FSI terms for 
finite Q^. Those FSI, induced by binary collisions between the struck and a core nucleon, 
depend on effective A^A^ interactions, matter properties, none of which are known with great 
precision. 

A first application is for inclusive cross and the agreement with data is excellent for NI 
parts and also for the immediate NE side x > 1 of the QEP. It deteriorates for the most 
elastic regions x > 1 for all angles where cross sections are smallest. We found that ternary 



collision contributions do not account for the above discrepancies. However, different single- 
nucleon momentum distributions easily cover the range of the problematic parts of the cross 
sections. 

We performed a scaling analysis for ratios ^^^'^^{XjQ"^) of cross sections for different 
nuclei and observed the predicted A-independence over the entire NE and NI regions. A 
second ratio i^'^ iyGiQ"^) was found to approach a plateau for all yc < with negative 
slope, decreasing for increasing Q^. We also studied FSI and concluded that those are only 
weakly g-dependent. Consequently the plateau contains the asymptotic limit as well as FSI. 
The same seems also to be compatible with results of a direct fit of scaling functions. Should 
the above be confirmed, there may be less incentive to study nuclear scaling in the high-Q^ 
regime as a tool to extract nucleon momentum distributions. From the agreement between 
predictions for d'^a and data (Fig. 1) we conclude that a weakly g-dependent FSI, which 
does not contradict any principle, is empirically correct. 

In our conclusion we assemble results obtained for a variety of inclusive scattering observ- 
ables, which cover a wide range of kinematic ranges, and which all are related to computed 
nuclear form factors. In addition to the above mentioned inclusive cross sections, we men- 
tion R ratios and ratios of moments of F^{x, Q^^. Generally there is agreement with data, 
which frequently span broad ranges of variation. All prediction use the underlying relation 
(§) between structure functions of a nucleus and of nucleons. Of course, the agreement 
does not prove the relation but neither does it disclose definite fiaws. In fact, one wonders 
whether (j^) is the result of some effective theory which, at least for the chosen observables, 
successfully replaces the fundamental QCD. 

The authors thank John Arrington for having communicated data and Brad Fillipone 
for discussions and correspondence on the subject matter. 
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Figure captions 

Fig. 1 Data0 and predictions^ for the CEBAF 89-001 experiment. 

Fig. 2. The NE part of the GRS-type scaling function ^^^'^^^^^yc < 0,Q^) © as 
function of Q^. Dots connect extracted values, drawn lines data, with calculated NI parts 
removed. 
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Table I 



In, \ ( r\T\ 


U 


X 


/02 ir^r^\T'i\ 
Q (GeV ) 


^C,Fe 


^Fe,Au 




23 


2.30 


2.26 


0.81 


1.03 


-0.4 


30 


1.95 


3.38 


0.70 


0.84 




45 


1.67 


5.46 


0.97 






15 


2.49 


1.05 


0.82 


1.00 


-0.2 


30 


1.37 


3.09 


0.98 


1.19 




55 


1.30 


5.78 


0.87 


1.24 




15 


1.02 


0.97 


1.18 


1.05 


0.0 


30 


1.01 


2.79 


1.04 


1.16 




74 


1.01 


5.77 


1.28 


0.84 




15 


0.65 


0.91 


0.97 


1.02 


0.2 


30 


0.72 


2.43 


1.00 


1.10 




74 


0.74 


4.54 


1.08 




0.4 


15 


0.43 


0.83 


1.00 


1.03 



Selection of cross section ratios Eq. (H). For each selected, narrowly-binned {yc), 

available data of the ratios are given for the smallest, some medium and largest (x, Q^) in 
the data sets. 
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